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INTRODUCTION : 


Advancements  on  technology  increase  the  potential  for  precision  in 
radar  navigation  and  communication  systems,  but  daily  variability  of  the 
ionosphere  continues  to  be  a  constraint  on  achievable  system  performance. 
Numerical  maps  which  provide  monthly  median  correction  for  ionospheric 
effects  have  been  derived  from  a  world-wide  climatology  of  ionospheric 
parameters;  their  use  alone  reduces  the  monthly  r.m.s.  residual  error  to 
about  20-25  percent  of  the  median  correction  in  daytime  and  30-35  percent 
at  night.  Now  20  percent  of  the  monthly  median  correction  for  ionospheric 
effects  is  often  still  the  dominant  error  component  in  systems  that  are 
free  to  calculate  performance  on  a  monthly  basis;  but  many  present  systems 
cannot  tolerate  the  individual  departures  of  50  to  200  percent  that  are 
averaged  away  in  these  monthly  statistics.  For  these  precision  systems 
adapted  predictions  have  been  proposed  which  would  use  various  field  meas¬ 
urements  to  make  real  time  corrections  for  ionospheric  effects. 

The  local  ionosphere  in  the  system  coverage  area  is  monitored  through 
measurements  of  ionospheric  parameters  such  as  the  critical  frequency  of 
the  F-region  (foF2),  total  electron  content  (TEC),  the  gradient  in  the 
lower  F-region;  or  through  measurements  derived  from  the  system  itself, 
such  as  the  signal  to  noise  ratio  measured  at  a  distant  site  and  retrans¬ 
mitted  to  the  system  control.  Some  system  users  may  not  have  the  capabili¬ 
ty  of  implementing  computer  models  based  on  world-wide  climatology,  and  in¬ 
stead  may  use  locally  monitored  TEC  or  foF2  to  estimate  monthly  median 
behavior  and  short  term  persistence  effects.  That  is,  they  produce  means 
of  several  days  to  predict  the  subsequent  hourly,  daily,  weekly  response 
of  their  system  to  the  current  state  of  the  ionosphere.  If  the  user  is 
only  concerned  with  ionospheric  behavior  in  the  area  being  monitored  by 
these  measurements,  such  a  technique  may  be  sufficient,  but  if  extrapola¬ 
tion  in  space  is  necessary,  the  use  of  a  model  is  required  to  account  for 
latitudinal  and  local  time  gradients. 

Adaptive  procedures  that  use  real  time  observations  to  reduce  average 
monthly  r.m.s.  errors  in  predictions  have  been  examined  using  foF2  data*, 
and  TEC  data^’3,4.  These  include  TEC  as  derived  from  measurements  of  Dop¬ 
pler  shift  on  ensuing  passes  of  the  Navy  Navigation  Satellites  (NNSS) ,  and  from 


1.  Wilkinson,  P.J.  (1970):  Prediction  Limits  for  foF2.  International 

Solar-Terrestrial  Predictions  Workshop  Proceedings.  April  1979, 
to  be  published  by  NOAA,  Boulder,  Colorado. 

2.  Donatelli ,  D.E.  and  R.S.  Allen  (1978):  Temporal  Variability  of  Iono¬ 

spheric  Refraction  Correction.  Proc.  IES  78,  Symposium  on  the  Ef- 
fect  of  the  Ionosphere  on  Space  and  Terrestrial  Systems,  Editor  J.M. 
Goodman,  NRL  and  ONR,  Washington,  D.C.,  1978:  490-496,  (/008-051-0069- 

1. 


Faraday  rotation  measurements  using  VHF  beacons  from  geostationary  satel¬ 
lites.  The  procedures  apply  to  any  type  of  synoptic  data,  and  rely  on  a 
reasonable  prediction  of  the  monthly  mean  of  the  observations,  the  standard 
deviation,  and  an  assumption  on  the  correlation  of  the  observations  in  the 
space-time  cell  being  considered.  Two  basic  procedures  will  be  discussed: 
the  first  addresses  multiplicative  processes  in  the  variability  of  the  ion¬ 
osphere;  the  second  addresses  additive  processes.  In  both  cases  the  results 
of  positive  and  negative  correlation  are  considered. 

An  evaluation  of  the  procedure  is  presented  considering  the  case  of 
multiplicative  processes.  A  baseline  mapping  of  monthly  mean  refraction 
correction  is  adapted  with  real  time  measurements  of  TEC;  the  procedure 
scales  the  map  of  a  system  coverage  area  with  a  ratio  of  observed  to  pre¬ 
dicted  TEC.  This  ratio  is  determined  at  a  reference  point  and  extended  in 
space  and  time.  Since  any  mapping  based  on  climatology  must  introduce 
model-dependent  bias,  an  ideal  baseline  mapping  is  approached  by  using  the 
actual  mean  of  the  observations  in  the  procedure. 

The  temporal  and  spatial  growth  of  the  r.m.s.  residual  error  is  exam¬ 
ined  using  the  concept  of  a  space-time  cell,  defined  as  the  extension  in 
space  and  time  from  a  reference  point  where  the  residual  error  is  defined 
as  zero,  to  the  limit  where  the  residual  error  of  the  adapted  map  is  equi¬ 
valent  to  the  variance  from  the  baseline  prediction.  For  an  individual 
system,  the  limit  would  be  modified  to  lie  within  the  error  budget  of  the 
system,  assuming  it  is  less  than  the  variance. 

Since  the  refraction  of  a  radio  signal  is  directly  proportional  to 
the  electron  content  along  the  ray  path,  a  system  user  may  translate  the 
results  to  be  presented  here  in  terms  of  system  parameters,  such  as  time 
delay  or  error  in  range  determination,  as  follows: 

Time  Delay  =  AT  =  40.3  x  TEC;  in  seconds 


40  3 

Range  Error  =  AR  =  cAT  =  — -4—  x  TEC;  in  meters 


3.  Leitinger,  R. ,  R.S.  Allen,  D.E.  Donatelli,  G.K.  Hartmann  (1978):  Adap¬ 

tive  Mapping  of  Ionospheric  Features.  Proc.  IES  78,  Symposium  on  the 
Effect  of  the  Ionosphere  on  Space  and  Terrestrial  Systems,  Editor  J.M. 
Goodman,  NRL  and  ONR,  Washington,  D.C.,  1978:  530-337,  #008-051-0069- 

1. 

4.  Donatelli,  D.E.  and  R.S.  Allen  (1979):  Ionospheric  Refractive  Correc¬ 

tion  using  an  Alaptive  Procedure.  International  Solar-Terrestrial 
Predictions  Wo^Kshop  Proceedings,  April  1979,  to  be  published  by  NOAA, 
Boulder,  Colorado. 
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where:  f  is 
c  is 
TEC  is 


the  frequency  of  the  radio  signal 
the  speed  of  light  in  a  vacuum 
in  units  of  electrons/meter ^  x  10^. 


ANALYSIS 


The  procedure  was  developed  from  the  assumption  that  the  major  compo¬ 
nents  of  the  daily  variability  of  the  ionosphere  could  be  described  as  large 
scale,  large  magnitude  fluctuations  that  are  slowly  varying,  and  therefore, 
reasonably  correlated  in  space  and  time.  These  fluctuations  result  from 
either  multiplicative  or  additive  processes;  the  analysis  is  similar  in 
each  case.  The  r.m.s.  error  when  using  the  adaptive  procedure  is  determined 
as  follows: 

£(°„  „  ~  PR  J 

k,  t  N  •  (1) 


where:  0R  t 

N 

t,R 


the  observation  at  the  point  where  the  correction  is  applied 

the  prediction  at  the  point  of  correction  using  an  observa¬ 
tion  from  a  reference  point 

the  number  of  observations 

the  time  and  place  of  the  applied  correction,  respectively. 


This  can  be  expanded  to: 

2  2  2  2 

XR,t  =  °R,t  +  MR>t  +  iPR>t  -  2  |-0RjtPR>t  (2) 

2 

Assuming  a  normal  distribution  with  mean  p  and  variance  a  ,  and  large  N  such 
that  N  '  N  -  1.  PR(t)  is  defined  by  the  process  being  considered,  as  fol¬ 
lows  : 


Case  1:  Multiplicative  process 

a.  positive  correlation 
0 


R,t 


,,  R,t 
vo,o 


negative  correlation 
0 


V  ’ 2 


*'r  ,  t 


Jo  ,0 
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Case  2:  Additive  process 


a.  positive  correlation 

PR,t  =  (°o,o  ~  '!o,o)  +  pR,t 

b.  negative  correlation 


PR,t  “  (mo,o  ~  °o ,  o  ^  +  Mu,t 


Vi  and  o  denote  the  mean  and  the  standard  deviation  of  the  observations,  re¬ 
spectively.  The  subscripts  0,0  and  R,t  denote  the  point,  in  space  and  time, 
of  reference  and  Qf  applied  correction,  respectively. 

To  evaluate  multiplicative  processes,  first  consider  case  la,  where 
by  substituting  for  PR  t  and  normalizing  to  the  moan,  equation  (2)  becomes: 


x;  <  _ 

R, t  _  R;t  +  o,o  _2 

„2  ,.2  72 


R,t 


R,  t 


Jo,o 


i.0  0 

o,o  R,  t 


Nyo,o!1R,t 


(3) 


From  the  definition  of  the  correlation  coeff iceient,  p  (t),  it  can  be  shown 
that:  R 


o,oMR,t 

Equation  (3)  can  then  be  written: 


SO  0  0  °R  t 

o»°  R, t  _  !  =  o,o  ’ 

Np 


VI  Ur,  t  R’t 
Mo,o  R,t 


(4) 


X  _  (J  U 

R,t  _  R.t  +  o,o  _  2  °»° 

uR,t 


0,0 


P  u  R.t 
o,o  R,  t 


(5a) 


Then  by  defining: 


o  =  q°»P  or  qo,o  =  a°R,t 

yo,o  |JR , t  ^0,0  UR,t 


and  substituting,  equation  5a  becomes: 


R.t 


JR,t  (1  +  °2  ~  2apR,t) 


(5b) 


If  a  negative  correlation  is  known  to  exist,  then  P  for  case  lb  is 
substituted  in  equation  (2)  with  the  following  result: 
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(6a) 
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T 


immtm 


Since  a  and  8  can  never  be  negative  then  |p  |  >  0.  If  p  =0, 
then  *R,t  >  °R>t;  if  «  =  3  =  0,  then  xR>t  =  oR>t.  *Jn  thl>se  cases’^e 

procedure  should  not  be  used. 

In  summary,  for  the  multiplicative  process,  the  procedure  is  effective 
if  U1  >  |P.  I  >  a/2  >  0 

K,  L 

and  for  the  additive  process,  if 

Ul  >  |pR>tl  >  3/2  >  0. 

Now  consider  the  ideal  cases  of  perfect  correlation,  pR  =  1  (case  a) 

and  anti-correlation,  p  =  -1  (case  b) :  ’ 

R,  t 

For  cases  la.  and  lb,  equations  5a  and  6a  reduce  to: 


*R,t  = 
pR,t 


o.o 


PR,t  ho,o 


(9a) 


and  5b  and  6b  reduce  to: 

,  =  o„  I 

R>c  R, t 

then  the  procedure  is  successful  within  the  limits: 
0<  ^£-<2  a*2± 


Jo,o 


JR,t 


0  <  a  <  2.- 

For  cases  2a  and  2b,  equations  7a  and  8a  reduce  to: 


(9b) 


XR  =  |o  ~  o'  I 

R,t  R, t  o,o 


(10a) 


and  7b  and  8b  reduce  to: 


XR,t*°R,t  11  ‘S| 


(10b) 


setting  the  limits  as: 


0  <  o  <  2o 

o,o  R,  t 


or  0  <  8  <  2. 
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EVALUATION 


From  this  analysis  it  can  be  seen  that,  if  perfect  correlation  or 
anti-correlation  exists,  the  success  of  the  adaptive  procedure  depends  on 
the  absolute  difference  in  the  percentage  standard  deviation  with  respect 
to  the  mean  between  the  reference  point  and  the  point  of  applied  correction 
("point"  referring  to  a  location  in  both  time  and  space) ,  for  a  multiplica¬ 
tive  process,  equation  9a. 

For  an  additive  process  it  is  simply  the  difference  in  the  standard 
deviation  between  points,  equation  10a.  Then,  it  is  possible,  by  comparing 
the  absolute  and  percentage  standard  deviations  between  relevant  points, 
to  estimate  the  potential  success  of  an  adaptive  procedure,  and  even  to 
determine  whether  .multiplicative  or  additive  processes  dominate.  For  this 
type  of  evaluation  to  be  meaningful,  there  should  be  simultaneous  data 
available  from  several  locations  within  a  typical  system  coverage  area, 
encompassing  diurnal,  seasonal  and  solar  cycle  variations. 

The  TEC  data  for  this  study  were  obtained  from  measurements  of  the 
Faraday  rotation  of  signals  from  beacons  of  geostationary  satellites^. 

The  available  data  were  in  two  sets:  the  first,  from  four  stations  for 
1968-69;  the  second,  from  three  stations  for  1974-76.  The  location  of 
monitoring  sites,  the  geostationary  satellite  that  was  monitored,  and  the 
location  of  the  420  km.  sub-ionospheric  point  along  the  ray  path,  are  pre¬ 
sented  in  Table  £. 

TABLE  I.  LOCATION  OF  TEC  MEASUREMENTS 

Station  •  Satellite  Sub-ionospheric  Point 


1968  -  69 

Hamilton,  Massachusetts  (HAM) 

ATS -3 

39N.  289E. 

Urbana,  Illinois  (URB) 

ATS-3 

36N.  274E. 

Stanford,  California  (STA) 

ATS-1 

34N.  234E. 

Edmonton,  Alberta  (EDM) 

1974  -  76 

ATS-1 

47N.  239E. 

Goose  Bay,  Labrador  (GSB) 

ATS-3 

48N.  289E. 

Hamilton,  Massachusetts  (HAM) 

ATS-3 

39N.  289E. 

Kennedy  Space  Flight  Center, 

Florida  (KSF) 

ATS-3 

26N.  280E. 

5.  Titlieridge,  J.E.  (1972):  Determination  of  Ionospheric  Electron  Content 
from  the  Faraday  Rotation  of  Geostationary  Satellite  Signals.  Planet¬ 
ary  and  Space  Science,  20:  353-369. 


The  first  set  of  data  represent  the  recent  solar  maximum,  and  the 
second,  the  subsequent  solar  minimum.  These  data  provide  a  good  basis  for 
evaluation  of  the  average  behavior  of  this  type  of  procedure  since  this 
solar  cycle  closely  approximates  the  average  of  the  past  20  solar  cycles 
with  Rz  =  110,  S  =  155  at  maximum,  and  Rz  =•  10,  S  =  71  at  minimum,  where 
Rz  and  S  are  the  twelve-month  running  mean  sunspot  number  and  solar  flux 
at  2800  Mhz,  respectively.  To  estimate  TEC  for  other  values  of  Rz  (or  S) 
it  is  still  possible  to  interpolate  or  extrapolate  to  a  reasonable  approxi¬ 
mation  by  assuming  a  linear  relationship  between  it.,  (or  S)  and  TEC.^ 

To  estimate  the  potential  of  adaptive  procedures,  first  consider 
multiplicative  processes,  as  defined  by  cases  la  and  lb  with  p  =  1  or  p  =  -1 
respectively.  In  this  case  the  percentage  r.m.s.  error,  Xr  t^R  t ’  esti- 
mated  by  equation  .9a,  which  requires  knowledge  of  the  percentage’standard 
deviation  for  the  reference  point  and  the  point  of  applied  correction.  This 
percentage  was  calculated  for  each  month  of  data  at  each  location  in  Table  I 
The  results  are  shown  in  Figures  la,b,c,d  for  the  months  of  January,  April, 
July  and  October,  respectively,  for  the  years  1969  and  1975,  allowing  exami¬ 
nation  of  diurnal,  seasonal,  and  solar  cycle  effects.  For  an  operational- 
type  assessment  it  seemed  more  appropriate  to  plot  the  curves  in  Universal 
Time.  From  these  curves  it  can  be  determined  between  which  locations  and 
during  which  times  of  day,  season  and  solar  cycle,  the  minimum  conditions 
are  met  for  using  the  adaptive  procedure. 

The  growth  of  the  percentage  error  in  time,  at  any  location,  may  be 
estimated  by  taking  progressive  differences  along  the  curve  for  a  particular 
location  from  some  choice  of  reference  time  on  that  curve.  When 
a  !  a. 


0, t  _  0,0  =  0,t 


the  limit  of  a  time  cell  is  defined.  In  the  same 


1  o.t  "0,01  1*0,  t 

manner,  the  limit  for  a  spacetime  cell  is  defined  by  taking  progressive 
differences  in  time  and  location  from  some  defined  reference  point,  to  the 


point  where 


o  o 

R,  t  _  0,0 


a 

=  R,t 


One  of  the  limitations  on  space-time  cell  extent  that  may  be  deter- 
minted  from  these  curves,  is  that  large  errors  can  be  generated  when  a 
scaling  factor  from  nightime  is  used  across  the  sunrise  terminator  to  adapt 
daytime  predictions.  This  is  a  greater  concern  at  solai  minimum  than  at 
solar  maximum  since  the  percentage  variability  fluctuates  more  radically; 
with  the  lower  mean  values  at  solar  minimum  a  small  absolute  variation  of 
TEC  will  produce  a  large  percentage  var  ation,  particular! v  at  nighttime. 


6.  DuLong,  D.D.  (19"’/);  Reduction  of  the  Uncertainty  of  Radar  Range 
Correction.  AFGL-TR-7 7-0125 ,  DDC#  ADA  046166. 


In  general,  the  extended  space-time  cells  that  are  possible  at  solar  maximum 
must  be  reduced  at  solar  minimum.  However,  the  remaining  variability,  in 
TEC  units  (electrons/meter2  x  10^) t  is  still  less  at  solar  minimum.  Note 
that  it  is  y^_,  not  the  percentage  y  /p  t  il"13!  of  primary  concern  to 
systems.  ’  ’ 

When  additive  processes  are  considered,  the  potential  of  the  adaptive 
technique  may  be  estimated  using  figures  2a,b,c,d.  These  are  the  curves 
for  the  standard  deviation  in  TEC  units  for  the  months  and  locations  of 
figure  1.  Again,  the  process  is  considered  for  p  =  1  or  p  =3  -1,  where  the 
error,  Xr  t  is  defined  by  equation  10a.  From  these  curves  Xr  t  i-S  °i>tain_ 
ed  directly  in  TEC  units  by  taking  differences  between  points, ’analogous  to 
figure  1.  This  figure  indicates  that  variability  in  terms  of  additive  pro¬ 
cesses  may  provide  a  better  description  of  nighttime  and  solar  minimum  con¬ 
ditions,  but  woul,d  be  less  descriptive  of  daytime,  particularly  at  solar 
maximum. 

.ADAPTIVE  PROCEDURE 


Since  the  daytime  <error  at  solar  maximum  causes  the  severest  system 
problems,  it  was  decided  that  this  study  should  examine  the  adaptive  pro¬ 
cedure  for  a  multiplicative  process  in  detail,  specifically  case  la.  The 
assumption  here  is  that  within  the  available  data  base,  the  correlation  is 
expected  to  monotonically  decrease  to  zero  with  distance,  and  negative  cor¬ 
relation  is  not  anticipated.  Based  on  correlation  studies^'®  this  is  a 
reasonable  assumption,  particularly  at  mid-latitudes. 

The  procedure  is  initiated  by  obtaining  a  measurement  of  TEC  at  a  ref¬ 
erence  point  (0Q  „)  determining  a  scaling  factor  from  the  ratio  of  the  ob¬ 
servation  to  the’baseline  prediction  for  the  reference  point  (Mq  q)  and 
scaling  the  baseline  map  over  the  coverage  area  for  a  succeeding’ time  period. 
The  adapted  prediction  (PR  t)  at  a  remote  time  and/or  site  is  then: 


P 


R,t 


0 


M 


0,0 

0,0 


M 


R,t 


where  MR  t  is  the  baseline  prediction  at  the  remote 


time  and/or  site.  As  noted  previously,  the  actual  monthly  mean  is  used  for 
the  baseline  prediction  at  each  site  to  remove  any  model-dependent  error 
from  the  ensuing  results.  Thus,  From  case  la: 


0 


Vt  -  -&1O 


||  •  The  TEC  dn^a  are  reduced  at  15-minute  intervals.  The 

0,0  ,<,C 


7.  Rush,  C.M.  (1076):  An  Ionospheric  Observation  Network  for  use  in  Short- 

Term  Propagation  Predictions.  Telecommunication  Journal,  43  (VIII): 
544-549.  . . 

8.  Klohuchar,  J.A.  and  J.M.  Johanson  (1977):  Correlation  Distance  of  Mean 

Daytime  Electron  Content.  AFGJ •-TR-77-0 L85 ,  DDC/I  ADA  048117. 
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adapted  prediction  is  applied  at  these  intervals  over  several  hours,  at 
all  locations  having  simultaneous  data.  Comparisons  are  made  with  the 
corresponding  observations,  and  the  r.m.s.  error,  xR  t  is  computed  each 
month  at  each  location,  at  each  15-minute  increment  from  the  same  reference 
times  each  day. 

TIME  CELL  RESULTS 

The  potential  of  this  simple  adaptive  technique  for  reducing  the 
r.m.s.  error  of  a  refractive  correction  is  demonstrated  in  Figure  3,  using 
data  from  Hamilton,  Massachusetts  for  January,  1969.  For  this  case  the 
local  ionospheric  measurements  are  used  to  predict  their  subsequent  tem¬ 
poral  behavior,  so  there  is  no  spatial  component  in  the  results.  In  figure 
3a  the  observed  10-day  mean  is  used  to  map  into  future  times;  in  figure  3b 
a  standard  AFGL  model  based  on  the  ionospheric  climatology  contained  in 
world  maps^  is  the  baseline.  Note  that  the  diurnal  shape  derived  from  the 
AFGL  model  differs  from  the  observed  ionospheric  behavior.  The  outer 
scale  on  the  figures,  to  be  used  for  the  means,  is  four  times  greater  than 
the  inner  scale,  which  is  to  be  used  for  the  standard  deviations  with  res¬ 
pect  to  the  means  and  for  the  r.m.s.  error  of  the  adapted  predictions  with 
respect  to  the  observations.  Although  there  is  up  to  20  percent  difference 
between  the  AFGL  predicted  and  the  observed  mean,  the  monthly  standard  de¬ 
viation  with  respect  to  each  mean  is  comparable.  Thus  on  the  basis  of  using 
a  monthly  prediction  without  adaption,  the  world-wide  monthly  median  clima¬ 
tology  predicts  as  well  as  the  mean  itself. 

The  growth  of  the  r.m.s.  error  of  the  adapted  prediction,  xR  t  is  de¬ 
scribed  by  the  thin  lines  originating  at  a  convenient  2  hour  internal,  where 
X..  n  In  3a‘,  xH  ^  extends  in  time  from  t  =  t  to  t  =  t  +  3  hours;  in 

o 

3b  from  t  =  0  to  t  =  tQ  and  held  constant  over  the  3  or  12  hour  time  cell; 
only  he  scaling  of  the  mean  (p  )  or  the  predicted  mean  (M  )  is  prog¬ 
ressed  with  time.  ">t 

The  dashed  lines  define  the  residual  error  for  a  3  hour  time  cell 

((x„,  t  +  3h)  and  indicate  that  a  scaling  factor  that  is  3  hours  old  is  use- 
H  o 

ful  throughout  the  daytime  hours,  unless  it  was  determined  prior  to  sunrise 
and  used  thereafter.  The  difference  between  la  and  lb  is  the  model-depend¬ 
ent  error,  most  evident  in  the  sunrise  area.  Comparing  these  curves  with 
the  standard  deviation  (oR  t)  with  respect  to  the  monthly  mean  and  the  stand¬ 
ard  deviation  (S  )  with  fespect  to  the  model  prediction,  it  is  apparent 
H,  t 


9.  Barghausen,  A.F.,  J.W.  Finney,  L.L.  Proctor,  L.D.  Schultz  (1969):  Pre¬ 
diction  Long-Term  Operational  Parameters  of  High-Frequency  Sky-Wave 
Telecommunications  Systems.  ESSA  Technical  Report  ERL  110-ITS-78. 
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that  the  time  cell  in  which  an  adapted  prediction,  based  on  a  single  scaling 
factor,  is  useful,  is  bound  by  the  reference  time  on  one  end  and  the  succeed¬ 
ing  solar  terminator  on  the  other.  (An  analysis  of  t  for  case  la,  that 
includes  the  effect  of  the  model  bias  can  be  found  in  Appendix  A. 

The  ana  Lysis  in  this  report  has  assumed  the  data  follows  a  normal 
distribution.  In  reality,  however,  it  is  somewhat  skewed  toward  higher 
values  in  daytime.  Therefore,  a  model  which  tends  to  predict  a  mean  that 
is  slightly  greater  than  the  observed  mean  in  daytime  will  often  succeed  in 
reducing  the  monthly  variance  at  that  time. 

In  a  previous  study^  the  performance  of  this  adaptive  technique,  pre¬ 
dicting  only  in  time  for  a  fixed  location,  was  examined  using  all  of  the 
archive  data  of  TEC  from  Hamilton,  Massachusetts,  1969  through  1976.  Table 
2  summarizes  these  results  for  the  daytime  and  nighttime  maxima,  averaged 
over  season,  at  both  solar  maximum  and  minimum,  in  TEC  units  of  electrons/ 
meter^  x  10^6.  This  table  shows  that  an  adapted  prediction  can  reduce  the 
residual  error  by  60  percent  in  a  one  hour  time  cell  and  still  by  as  much 
as  30  percent  in  a  three  hour  time  cell  for  daytime  solar  maximum. 

The  time  cell  for  useful  adapted  predictions  is  determined  by  the  cor¬ 
relation  times  of  fluctuations  in  TEC.  Since  the  basic  data  for  this  study 
were  reduced  at  15  minute  intervals,  fluctuations  with  periods  less  than  30 
minutes  are  not  observable.  Their  amplitude  is  assumed  to  be  much  smaller 
than  the  longer  period  variations  that  are  being  tracked  here.  The  present 
simple  adaptive  technique,  predicting  persistence  of  a  direction  of  varia¬ 
tion  for  times  of  15  minutes  to  three  hours,  tends  to  mitigate  against  just 
those  long  period  TF.C  variations  which  seem  to  be  a  major  component  at  solar 
maximum.  Thus  this  adaptive  procedure  is  most  successful  at  solar  maximum, 
yet  the  residual  is  of  the  same  absolute  magnitude  as ? the  residual  at  solar 
minimum.  This  suggests  that  the  short  period  variations  of  TEC  predominate 
at  solar  minimum,  and  are  still  present  with  the  same  amplitude  at  solar 
maximum,  masked  by  the  characteristic  larger  scale  variations. 

During  sunrise  and  sunset  periods,  the  effects  of  production  and  loss 
create  gradients  that  dominate  all  other  fluctuations.  Similar  gradients 
may  occur  during  magnetic  disturbance,  particularly  in  the  region  of  the 
auroral  zone  and  the  trough  in  high  and  mid-latitutdes.  The  rate  at  which 
these  changes  occur  would  determine  the  time  cells  for  these  periods.  This 
is  illustrated  in  figure  2  where  the  observed  TEC  is  compared  to  the  results 
of  using  adapted  predictions  on  a  highly  disturbed  day  at  Hamilton,  Mas¬ 
sachusetts. 

Figure  4a  compares  the  prediction  for  the  month,  y  the  observed 
TEC  for  the  day,  0„  the  adapted  prediction  for  a  30  min&te  time  cell, 

p  " ,  t 

H,t  +  30m  and  a  3-hour  time  cell,  P„  .  Figure  4b  shows  the  differ- 
°  H,  >:  +  Jh. 


6.  DuLong,  D.D.  (1977):  Reduction  of  the  Uncertainty  of  Radar  Range 
Correction.  AFGL-TR-77-0125 ,  DDC//  ADA  046166. 
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TABLE  2.  RESIDUAL  ERRORS  IN  TEC  PREDICTION  TIME  CELLS 


WINTER  VERNAL  EQUINOX  SUMMER  AUTUMNAL  EQUINOX 


DAY 

NIGHT 

DAY  NIGHT 

DAY 

NIGHT 

DAY 

NIGHT 

Solar 

Maximum 

<« z  = 

110;  S 

=  155) 

“H.t 

34.4 

5.7 

44.8 

8.2 

25.4 

8.2 

41.1 

6.7 

°H,t 

6.0 

1.8 

7.8 

2.5 

4.0 

2.2 

7.8 

2.2 

XH,t0+3h 

4.0 

1.5 

4.5 

1.5 

2.7 

1.5 

4.0 

1.8 

XH,to+lh 

2.2 

1.0 

2.5 

1.0 

1.6 

1.2 

1.6 

0.7 

XH,to+30m 

1.5 

0.6 

1.5 

0.6 

1.2 

0.6 

1.3 

0.6 

Solar 

Minimum 

<*Z  = 

10;  S  = 

71) 

PH,t 

10.5 

1.5 

10.5 

1.5 

9.0 

2.1 

11.2 

2.1 

°H,t 

1.8 

0.6 

2.5 

0.6 

1.5 

0.5 

1.5 

0.6 

XH, t  +3h 
o 

1.8 

0.9 

2.1 

0.5 

1.5 

0.5 

2.1 

0.6 

XH,tc+lh 

1.2 

0.3 

1.2 

0.1 

0.6 

0.3 

1.2 

0.3 

XH,tQ+30m 

0.9 

0.1 

0.6 

0.1 

0.4 

0.1 

0.6 

0.1 

The  symbols  used  represent  the  following: 

u  -  The  mean  value  of  TEC 

H,  t 

Ojj  t  -  The  r.m.s.  deviation  of  TEC  from  the  mean,  representing 

the  day-to-day  variability  of  the  ionosphere 

V 

H,tQ+3h  -  The  residual  error  using  an  adapted  prediction  in  a 
3-hour  time  cell 


XH,t  +lh  "  The  residual  error  using  an  adapted  prediction  in  a 
1-nour  time  cell 

X„  ,-n  -  The  residual  error  using  an  adapted  prediction  in  a 

’ Co  m  30-minute  time  cell. 
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ences  between  the  observed  TEC  and  the  predictions  in  each  case.  The  error 
using  the  3-hour  time  cell  0„  .  -Pu  .  , I  is  equivalent  in  magnitude 

n  j  t  oh  n  ,  t  qi  j  n  | 

to  that  using  the  monthly  mean  |oR  t|,but  large  excursions  have  been 

introduced  at  time  periods  when  they  would  otherwise  not  exist.  The  adapted 
prediction  using  a  30  minute  time  cell,XR  +30m  Prova®es  a  60  percent  re¬ 
duction  in  the  maximum  error  despite  the  steep  gradients  in  TEC  on  this  dis¬ 
turbed  day. 

SPACE-TIME  CELL  RESULTS 

The  preceding  results  assess  the  adapted  prediction  when  applied  at 
the  reference  location,  progressed  in  time,  thus  including  only  a  temporal 
component  of  variability.  A  more  valid  assessment  of  this  procedure  requires 
inclusion  of  the  spatial  component  as  well.  Spatial  variations  have  been 
previously  examined®* 10* H  using  data  from  the  TRANSIT  satellites  reduced  to 
provide  TEC  along  the  satellite  path,  and  TEC  measurements  from  several 
closely  spaced  monitors  of  the  ATS-6  geostationary  satellite.  The  variation 
of  ionospheric  features  and  some  of  the  implications  of  attempting  to  track 
these  features  in  space  and  time  were  considered  and  will  be  further  examined 
here. 

The  following  figures  represent  the  actual  impact  on  systems  of  the  use 

of  this  procedure.  It  will  be  seen  that,  qualitatively  the  results  agree 

with  the  estimates  inferred  in  figure  1.  Figures  5-8  present  the  1969  data 

for  the  months  of  January,  April,  July  and  October,  respectively,  and  the 

results  of  its  use  in  the  adaptive  procedure;  figures  9-12  present  the  1975 

data  for  the  same  months  in  the  same  manner.  The  (a)  figures,  5-12,  compare 

the  mean,  p  standard  deviation, o  and  the  residual  error 

R, t ,  R, t , 

minute,  (XR> to+15m> ,  1-hour  (xR>ty+lh),  and  3-hour  (X^^)  time  cell  at 

each  of  the  stations  in  each  set  of  data.  The  (b)  figures,  5-12,  compare 
the  standard  deviation  (o  )  for  each  station  with  the  residual  error  for 

a  15-minute  (xR^+15m)  1-hour  (xR>t  +lh>  and  3-hour  (Xr^^)  space-time 

cell,  the  spatial  reference  being  one  or  two  of  the  other  stations  within 
the  set  for  the  year. 

8.  Klobuchar,  J.A.  and  J.M.  Jonanson  (1977):  Correlation  Distance  of  Mean 
Daytime  Electron  Content.  AFGL-TR-77-0185,  DDC//  ADA  048117. 

10.  Allen,  R.S.  (1977):  Considerations  Relative  to  Adapting  TRANSIT  Observa¬ 

tions  to  Predicting  Radar  Range  Corrections.  AFGL-TR- 7 7-0004 ,  DDC# 

ADA  038238. 

11.  Davies,  K. ,  Degenliardt,  G.K.  Hartmann,  R.  Leitinger,  (1977):  Electron 

Content  Measurements  over  the  U.S.,  Joint  Radio  Beacon  Program  NOAA/ 
MPAE/C.RAZ,  Station  Report  ATS-6°,  94°W. ,  Edited  by  Max  Planck  Institut 
fUr  Aeronomie,  Lindau,  West  Germany. 
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Since  the  local  time  variation  at  each  location  is  inherent  In  the 
mean,  the  variations  being  tracked  by  this  procedure  need  not  be  local  time 
dependent.  It  can  be  seen  that  the  results  using  a  15-minute  time  cell  when 
projecting  in  space,  are  similar  or  slightly  better  than  using  a  1-hour,  and 
consistently  better  than  using  a  3-hour  time  cell.  This  is  independent  of 
the  distances  and  directions  considered  here.  It  should  be  noted,  (see  fig¬ 
ure  3)  that  the  error  growth  rate  is  greatest  during  the  first  hour  and 
tends  to  decrease  thereafter.  The  actual  error  maximizes  and  minimizes 
where  the  changes  in  variance  are  zero,  near  the  midday  peak  and  nighttime 
minimum,  respectively. 

Direct  quantitative  comparisons  are  not  possible  between  the  1969  and 
1975  data  since  the  spatial  coherence  in  latitude  is  less  than  in  longitude. 
The  station-pairs  GSB-HAM  and  HAM-KSF  of  the  1975  data  have  latitudi¬ 
nal  separations  tha’t  are  equal  to,  or  greater  than,  their  longitudinal  sep¬ 
arations.  This  is  also  true  for  the  1969  STA-KDM  station-pair,  but  for 
the  HAM-URB  and  URB-STA  station-pairs,  the  longitudinal  exceeds  the  latitudi 
nal  separations.  Therefore,  only  the  1969  STA-EDM  station-pair  results  may 
be  qualitatively  compared  to  the  1975  results  for  an  evaluation  of  solar 
cycle  variation.  However,  the  1969  and  1975  Hamilton  results  may  be  quanti¬ 
tatively  compared. 

At  solar  maximum  an  improvement  of  50  percent  or  more  is  achievable  in 
daytime  within  most  of  the  space-time  cells  presented  here;  at  solar  minimum 
this  is  only  achievable  in  daytime  within  a  space-time  cell  with  a  1-hour 
temporal  limit.  Some  degree  of  improvement  is  also  possible  at  nighttime, 
solar  maximum;  but  at  solar  minimum  nighttime,  new  errors  may  be  introduced 
with  this  procedure.  This  was  anticipated  since  it  is  assumed  that  additive 
rather  than  multiplicative,  processes  dominate  at  nighttime  and  at  solar 
minimum. 

ESTABLISHING  LIMITS  FOR  A  SPACE-TIME  CELL 

An  approximation  for  space-time  cells  is  suggested  by  the  results  of 
correlation  studies^’®.  For  northern  mid-latitude  regions,  these  indicate 
a  ratio  of  1:2  between  latitude  and  longitude  correlation  distances  in  de¬ 
grees.  Therefore,  an  ellipse  with  an  axial  ration  of  1:2  could  be  used  to 
describe  the  spatial  extent  of  an  applied  correction.  While  the  axial  ratio 
would  remain  constant,  the  spatial  extent,  in  kilometers,  would  decay  in 
time.  The  initial  distance  would  be  set  by  the  desired  accuracy,  and  assum¬ 
ing  space-time  equivalence  in  longitude. 


7.  Rush,  C.M.  (1976)'  An  Ionospheric  Observation  Network  for  use  in  Short- 

Term  Propagatior  Predictions.  Telecommunication  Journal,  43  (VIII): 
544-549. 

8.  Klobuchar,  J.A.  and  J.M.  Johanson  (1977):  Correlation  Distance  of  Mean 

Daytime  Electron  Content.  AFGL-TR-77-0185 ,  DDC //  ADA  048117. 
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This  concept  may  be  valid  for  tracking  features  even  when  their  growth 
and  damping  rates,  in  space  and  time,  are  not  explicitly  defined.  The  fea¬ 
tures  observed  in  Figure  A,  which  are  apt  to  be  confined  to  a  particular 

local  time  and  magnetic  field  shell  sector,  are  an  example.  Note  the  error 
of  the  30-minute  time  cell  prediction  in  Figure  A.  If  this  error  is  within 
the  tolerance  of  a  hypothetical  system,  the  instantaneous  boundary  of  a 
space-time  cell  can  be  defined  by  a  30-minute  equivalent  longitude,  7.5 
degrees.  At  Hamilton  latitudes  this  would  represent  a  6A5  km  East-West 
axis  and  a  A20  km  North-South  axis.  If  a  correction  is  applied  at  a  later 
time,  this  ellipse  is  reduced;  for  example,  15-minutes  after  the  initial 
observation  the  East-West  axis  would  be  reduced  to  323  km  and  the  North-South 
to  210  km.  At  the  end  of  30-minutes  the  correction  would  only  apply  at  the 
location  of  the  initial  observation. 

The  general  case  is  seen  in  figures  5-8  where  a  reduction  in  variance 

of  50  percent  at  solar  maximum  in  daytime  can  be  achieved  within  a  3-hour 

time  cell,  or  equivalently,  A5°  longitude.  This  is  equivalent  to  the  URB- 
STA  station-pair.  An  observation  at  Urbana  can  be  used  instantaneously  to 
provide  an  adapted  prediction  for  Stanford,  but  this  spatial  extent  would 
decay  such  that  after  3  hours,  the  adapted  prediction  would  apply  only  to 
Urbana  itself. 

The  extent  of  the  space-time  cells  is  greatest  at  solar  maximum  day¬ 
time  when  large-scale,  long  period  variations  predominate  and  systems  are 
most  concerned  in  an  r.m.s.  sense.  For  systems  that  are  vulnerable  to 
individual  features,  such  as  those  of  figure  A,  the  extent  of  the  space- 
time  cell  must  be  modified  accordingly. 

CONCLUSION 

When  using  adapted  predictions  analysis  of  r.m.s.  error  must  be  based 
on  whether  the  dominant  features  are  created  by  additive  or  multiplicative 
processes,  and  whether  their  extended  effects  are  positively  or  negatively 
correlated.  From  this  analysis  it  was  determined  that  the  r.m.s.  error  for 
additive  processes  is  dependent  on  the  absolute  difference  in  the  standard 
deviation  of  the  two  space-time  points  under  consideration.  For  a  multi¬ 
plicative  process  the  absolute  difference  in  percentage  standard  deviation 
with  respect  to  the  mean  is  the  related  parameter.  A  comparison  of  these 
parameters  for  the  available  data  indicated  that  for  the  solar  maximum  per¬ 
iod,  when  ionospheric  refraction  effects  on  systems  are  the  most  severe, 
multiplicative  processes  dominate. 

Evaluation  of  the  adaptive  procedure  required  that  data  be  confined  to 
a  region  typical  of  a  system  coverage  area.  The  available  data  were  limited 
to  the  northern  mid-latitude  American  sector.  From  correlation  studies  it 
was  deemed  reasonable  to  assume  a  positive  correlation  throughout  the  region. 
Thus  it  was  decided  to  evaluate,  in  detail,  an  adaptive  procedure  that  con¬ 
sidered  multiplicative  processes  with  positive  correlation.  The  required 

condition  was  shown  to  be  1  >  p  >  n  >  0 

R,t  2 

where:  p  Is  the  correlation  coefficient  between  points  and 
K,  t 


a  is  the  ratio  of  the  percentage  r.m.s.  deviation  at  the  reference  point 
to  that  at  the  point  of  applied  correction. 

In  daytime  an  improvement  of  50  percent  is  readily  achievable:  at 
solar  maximum  within  a  3-hour  time  cell;  at  solar  simimum  within  a  1-hour 
time  cell.  At  nighttime,  the  percentage  r.m.s.  deviation  is  usually  greater 
than  daytime,  but  the  absolute  deviation  is  less;  therefore,  the  need  for 
correction  is  not  as  crucial.  However,  if  adapted  predictions  are  used,  a 
scaling  factor  obtained  at  nighttime  should  not  be  projected  across  the  sun¬ 
rise  terminator,  as  the  rapid  change  in  percentage  r.m.s.  deviation  during 
the  sunrise  growth  period  introduces  errors.  The  sunset  decay  period  is 
less  critical  since  the  percentage  deviation  usually  increases  at  nighttime, 
therefore  a  would  decrease,  allowing  a  greater  range  for  ,  thereby  in¬ 
creasing  the  probability  that  the  requirement  pR  t  >  n  met 

A  space-time  cell  has  been  defined  in  terms  of  an  ellipse  of  constant 
axial  ratio,  with  a  spatial  extent  that  decays  with  time.  • It  was  found  that 
space-time  cells  with  an  extent  of  several  thousand  kilometers,  or  several 
hours,  or  an  equivalent  space-time  combination,  are  operationally  useful  at 
solar  maximum.  At  solar  minimum  the  absolute  value  of  the  r.m.s.  deviation 
is  low,  therefore  the  need  for  correction  is  less  critical .  For  example, 
the  residual  error,  using  a  one-hour  time  cell  at  solar  maximum,  is  equi¬ 
valent  to  the  average  monthly  variance  at  solar  minimum.  Therefore,  the 
adaptiveprocedure  succeeds  in  reducing  the  error  at  solar  maximum  to  the 
lower  impact  solar  minimum  levels. 
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Appendix  A 


For  actual  system  use,  the  monthly  mean  ionospheric  behavior  would  be 
predicted  by  a  numerical  model  derived  from  world-wide  climatology.  Any 
model  has  some  inherent  bias  which  may  have  little  consequence  on  r.m.s. 
error  when  the  primary  concern  is  predicting  the  mean  behavior,  but,  as 
shown  in  Figure  3b,  this  bias  can  effect  the  adaptive  procedure.  The  fol¬ 
lowing  derivation  will  provide  some  perspective  on  the  import  of  this  bias. 


The  r.m.s.  error  is  defined  as: 


X 


R,t 


y 

N 


(Al) 


where  PR  t  is  the  prediction  of  the  adaptive  procedure.  For  the  case  of  a 

multiplicative  process  within  an  area  of  positive  correlation,  which  is  the 
only  case  we  will  examine  here: 


M  nR,t 
0,0 


(A2) 


where:  M  is  the  mean  as  predicted  by  the  model  and  the  subscripts  0,0  and 
R,t  denote  the  reference  point  and  the  point  of  applied  correction,  respec¬ 
tively,  in  space  and  time. 


Substituting  this  in  equation  Al  and  expanding: 

i2 


T.  0 


%t 


_  M2  >:02  M  EC)  0 

R,  t  k, t  o,o  _  2  R»  t  R, t  o,o 

N  2  N  ~  M  N 

M  o,o 

o,o 


(A3) 


then  substitute  o„  _  a  ,  and  p„  in  equation  (A3)  to  arrive  at: 
R,t,  o,o  R, t 


XR,t  =  °R,t  +  PR,t  +  -/’-t  (0?  +  p2  )  -  2  \t  ,  , 

M2  °’°  °’°  (aR,t°o,oPR,t  +  Po,o,JR,t 

0,0  0,0 

(A4) 


To  parallel  case  la  in  the  text,  normalize  to  the  observed  mean: 

.2 


*L±  =  °Lt  +  l  +  \t  ‘ojo  Ao.o  + 

l,R,t  UR,t  ujj.  M2  L? 

R, t  o,o  \  o,o 


A  -  2  ^QjO  r_ R,t  gq,p 


V 


R* t  Mo,o  \"R,t  Ho,o 


p  4  11 

R,  t 


(A5) 
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Figure  Captions 


Figure  ia  -  The  diurnal  curves  of  percentage  standard  deviation  in  TEC  for 

the  month  of  January  at  four  stations  in  1969  and  three  stations 
in  1975. 

Figure  lb  -  Same  as  la  but  for  the  month  of  April. 

Figure  lc  -  Same  as  la  but  for  the  month  of  July. 

Figure  Id  -  Same  as  la  but  for  the  month  of  October. 

Figure  2a  -  The  diurnal  curves  of  standard  deviation  in  TEC  for  the  month 

of  January  at  four  stations  in  1969  and  three  stations  in  1975. 

Figure  2b  -  Same  as  2a  but  for  the  month  of  April. 

Figure  2c.  -  Same  as  2a  but  for  the  month  of  July. 

Figure  2d  -  Same  as  2a  but  for  the  month  of  October. 

Figure  3a  -  Comparison  in  TEC  for  January  1969  at  Hamilton,  Massachusetts 
of  the  observed  mean,  standard  deviation,  r.tn.s.  error  of  a  3- 
hotir  time  cell,  and  error  growth  in  time  using  an  adapted  pre¬ 
diction  up  to  3  hours  after  a  reference  time. 

Figure  3b  -  Similar  to  3a  but  using  a  model  prediction  of  mean  TEC.  The 

error  growth  in  time  using  the  adapted  prediction  is  shown  for 
12  hours  after  the  reference  time. 

Figure  4a  -  TEC  on  a  magnetically  disturbed  day,  compared  to  the  observed 
TEC  for  the  day  with  the  mean  and  the  adapted  predictions  for 

the  day  using  a  30-minute  and  a  3-hour  time  cell. 

Figure  4b  -  The  residual  error  in  TEC  using  the  predictions  of  4a  for  the 
day . 

Figure  5a  -  The  diurnal  curves  in  TEC  for  January  1969  from  four  stations 
comparing  the  mean  standard  deviation  and  the  residual  error 
using  the  adapted  predictions  with  15-m.inute,  1-hour  and  3- 
hour  time  cells. 

Figure  5b  -  Comparing  the  residual  error  using  the  indicated  space-time 

cells  to  the  standard  deviation  of  the  stations  of  figure  5a. 

The  arrows  indicate  the  direction  of  reference  to  applied 
correction. 

Figure  6a  -  Same  as  5a  but  for  April  1969. 
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Figure  6b  -  Same  as  5b  but  for  April  1969 

Figure  7a  -  Same  as  5a  but  for  July  1969. 

Figure  7b  -  Same  as  5b  but  for  July  1969. 

Figure  8a  -  Same  as  5a  but  for  October  1969. 

Figure  8b  -  Same  as  5b  but  for  October  1969. 

Figure  9a  -  Same  as  5a  but  for  three  stations  in  January  1975 

Figure  9b  -  Same  as  5b  but  for  three  stations  in  January  1975 

Figure  10a  -  Same  as  9a  but  for  April  1975. 

Figure  10b  -  Same  as  9b  but  for  April  1975. 

Figure  11a  -  Same  as  9a  but  for  July  1975. 

Figure  lib  -  Same  as  9b  but  for  July  1975. 

Figure  12a  -  Same  as  9a  but  for  October  1975. 

Figure  12b  -  Same  as  9b  but  for  October  1975. 
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